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Abstract. We suggest an algorithm for derivation of the Picard-Puchs system 
of Pfaffian equations for Abehan integrals corresponding to semiquasihomoge- 
neous Hamiltonians. It is based on an effective decomposition of polynomial 
forms in the Brieskorn lattice. The construction allows for an explicit upper 
bound on the norms of the polynomial coefficients, an important ingredient in 
studying zeros of these integrals. 



1. Introduction 

Given a polynomial in two variables / G M.[x, y] and a polynomial 1-form lo on M^, 
how many isolated ovals S on the level curves / = const may satisfy the condition 



JS 



Lj = 07 This is the long-standing infinitesimal Hilbert problem, see [Arn94]. The 
answer is to be given in terms of the degrees of / and lu. 

A recent approach to this problem, suggested in [NY99, NYOl, YakOl] is based on 
the fact that periods of polynomial 1-forms restricted on level curves of polynomials, 
satisfy a system of differential equations with rational coefhcients, called the Picard- 
Fuchs system. Under certain restrictions on the monodromy group, the number of 
zeros of solutions of such systems can be estimated from above in terms of the 
magnitude of coefhcients of this system, more precisely, the norms of its matrix 
residues. Thus it becomes important to derive the Picard-Fuchs system for Abelian 
integrals so explicitly as to allow for the required estimates for the residues. 
In [NYOl] a Fuchsian system was derived in the hypergeometric form 

it-i + A)i = Bi, i^iHt), (1.1) 

where I(t) — {Ii{t), ■ ■ ■ ,Ii{t)) is a collection of integrals of some monomial forms 
over any oval of the level curve {/ = t}, and A, B are two constant {I x Z)-matrices of 
explicitly bounded norms, depending on / (1 always stands for the identity matrix 
of the appropriate size). The rational matrix function R{t) — {t ■ 1 + A)^^B has 
only simple poles and the norm of its matrix residues can be explicitly majorized 
provided that the eigenvalues of A remain well apart. This allows to solve the 
infinitesimal Hilbert problem for all polynomials / whose critical values (after a 
suitable normalization) are sufRciently distant from each other. What remains is 
to study the case of confluent critical values (including those at infinity). 

In a general hypergeometric system (1.1), the residues may or may not blow 
up as some of the singular points tend to each other. The particular feature of 



Date: January 2002. 

1991 Mathematics Subject Classification. Primary 34C08; Secondary 34M50, 32S20, 32S40. 
The research was supported by the Israeli Science Foundation grant no. 18-00/1. 



1 



2 



SERGEI YAKOVENKO 



the Picard-Fuchs system is its isomonodromy: the monodromy group remains the 
same under deformations of / (at least for sufficiently generic /). This implies that 
even if the explosion of residues occurs, it cannot be caused by the explosion of 
the eigenvalues. In order to find out what indeed happens with the residues, the 
first step is to write down as explicitly as possible the Picard-Fuchs system as a 
flat meromorphic connexion with singularities in the holomorphic bundle over the 
variety of all polynomials / of a given degree. 

This problem is solved in the paper for polynomials with a fixed principal 
(quasi) homogeneous part having an isolated critical point at the origin. 

As an auxiliary first step, we need to describe explicitly the structure of the 
relative cohomology module. While the subject is fairly classic and sufficiently 
well understood, the existing proofs do not allow for the quantitative analysis. We 
suggest an alternative, completely elementary construction that immediately yields 
all necessary bounds. This construction, exposed in §2 is based on "division by 
/", a lemma distilled from the paper [Fra88] by J.-P. Frangoise. The Pfaffian form 
of the Picard-Fuchs system is derived in §4. In the last section we mention some 
simple properties of the derived system and formulate a conjecture that it has only 
logarithmic singularities in the affine part. 

2. Relative cohomology revisited 

2.1. Relative cohomology, Brieskorn and Petrov modules. Denote by A'', 

fc = 0, 1, . . . ,n, the module of polynomial fc- forms on the complex affine space C" 
for a fixed n > 1. If / G C[xi, . . . ,a;„] ~ A" is a polynomial, then the collection 
df A A''~^ of fc-forms divisible by d/ e A^, is a C-linear subspace in A'', and the 
quotient 

A'} = A''/df AA''-\ fc = l,...,n, (2.1) 

is called the space of relative k-forms. Since the exterior derivative d preserves 
divisibility by df, the relative de Rham complex A*, 

— >A)-^A) U A^-i A^ 0, (2.2) 

naturally appears. A form u) € A*^ is called relatively closed \i duj = df /\r] and 
relatively exact if a; = df A^ + d9 for appropriate rj G A'' and ^,0 G A''~^. The rela- 
tive cohomology groups ffj = _ff'^(A*), relatively closed fc- forms modulo relatively 
exact ones, are important characteristics of the polynomial /. 

Together with the natural C-linear structure, the relative cohomology groups 
i?j possess the structure of a module over the ring C[/] = /*C[a;i, . . . ,a;„]. This 
follows from the identity 

f -{df AT] + d0)=df A ifr, -9) + d{f0). (2.3) 

meaning that relatively exact forms are preserved by multiplication by /. 

As is well-known, the highest module H^j, as well as all H'j with < fc < n — 1, 
is zero. Instead, we consider another important module, called Brieskorn module 
(lattice) [Bri70, DSOl, DouOl], defined as the quotient 



Bf = A"-/df AdA 



(2.4) 
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and the C[/]-module P/, the quotient of all (n — l)-forms by the closed (n — 1)- 
forms, 

Pf = A^'-^/idf A + D H]-\ (2.5) 

The latter is an extension of H^~^: the quotient Pf/H^~^ is naturally isomorphic 
to the finite-dimensional C-space = A"/df A In several sources, P/ is 

referred to as the Petrov module. The exterior differential naturally projects as a 
bijective map d: P f ^ B f which obviously is not a C[/]-modulc homomorphism. 

Clearly, a relatively exact (closed) form is exact (resp., closed) after being re- 
stricted on any nonsingular level set f~^{t) C C", t e C, since df vanishes on all 
such sets. 

The inverse inclusion is considerably more delicate. Gavrilov studied the case 
n = 2 and proved that for a 1-form with exact restrictions on all level curves 
f~^{t) C to be relatively exact, it is sufficient to require that the polyno- 
mial / has only isolated singularities and all level curves f~^{t) be connected 
[Gav98, Gav99]. This result generalizes the earlier theorem by Ilyashenko [Ily69]. 
A multidimensional generalization in the same spirit was obtained by I. Pushkar' 
[Pus97]. The affirmative answer depends on the topology of a generic level set 
f~^{t) (its connectedness for n = 2 or vanishing of the Betti numbers bk for k 
between and n - 2, see [DS93, BDOO]). 

Both the isolatedness and connectedness assumptions can be derived from a sin- 
gle assumption that the principal (quasi)homogeneous part / of the polynomial / 
has an isolated critical point at the origin: such polynomials are called semiquasi- 
homogeneous [AGV85]. For two variables with equal weights it suffices to require 
that / factors as a product of pairwise different linear homogeneous terms. 

2.2. Computation of relative cohomology. Besides the above question on the 
relationship between the algebraically defined cohomology of the relative de Rham 
complex and analytically defined cohomology of (generic) fibers, the natural prob- 
lem of computing H* arises. 

This problem was addressed in the papers [BDOO, DSOl, DS93, DouOl, Gav99, 
Gav98] mentioned above. Using analytic tools or theory of perverse sheaves and TP- 
modules, they prove that under certain genericity-type assumptions on /, the high- 
est relative cohomology module H'^~^ and the Petrov module P/ are finitely gen- 
erated over the ring C[/]. For semiquasihomogeneous polynomials one can describe 
explicitly the collection of generators for B / , the polynomial forms wi , . . . , € 
A^~^ such that any other form uj e A"'~^ can be represented as 
I 

io = Y^Pi^i + dfATj + d^, Pi= Mf) e C[/], 7?,^ G (2.6) 

i=l 

with appropriate polynomial coefficients pi that are uniquely defined. 

The proofs of this and related results, obtained in either analytic or algebraic 
way, are sufficiently involved. In particular, it is very difficult if possible at all to 
get an information on (i) how the decomposition (2.6) depends on parameters, in 
particular, if / itself depends on parameters, and (ii) how to place explicit quanti- 
tative bounds on the coefficients Pi{f) in terms of the magnitude of coefficients of 
the form lu. For example, to extract such bounds from the more transparent ana- 
lytic proof by Gavrilov, one should place a lower bound on the determinant of the 
period matrix of the forms over a system of vanishing cycles on the level curves 
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f~^{t). The mere nonvanishing of this determinant is a deUcate assertion whose 
proof in [Gav98] is incomplete (a simple elementary proof was supplied by Novikov 
[NovOl]). The explicit computation of this determinant for a specific choice of the 
generators tOi was achieved by A. Glutsuk [GluOO], but the answer is given by a 
very cumbersome expression. 

In the next section we suggest an elementary derivation of the formula (2.6) under 
the assumption that the polynomial / is semiquasihomogeneous. This derivation: 

1 . gives an independent elementary demonstration of the Gavrilov-Bonnet-Dimca 
theorem for the most important particular case of semiquasihomogeneous 

polynomials; 

2. proves that the polynomial coefficients pi and the forms rj, 9 from the decom- 
position (2.6) depend polynomially on the coefficients of the non-principal 
part of /, provided that the principal quasihomogeneous part of / remains 
fixed; 

3. yields the collection of the coefiicients (pi, . . . ,pi) of (2.6) as a result of appli- 
cation of a certain linear operator to the form ui. The norm of this operator 
can be explicitly bounded in terms of / (and the chosen set of generators 
{u)i}) and the degree degw. 

3. Bounded decomposition in the Brieskorn and Petrov modules 

3.1. Degrees, weights, norms. In this section we first consider quasihomoge- 
neous polynomials from the ring C[x] = C[a;i, . . . , with rational positive weights 
Wi = deg Xi normalized by the condition wi + ■ ■ ■ + Wn = n to simplify the treat- 
ment of the most important symmetric case when Wi = 1. The symbol deg / always 
means the quasihomogeneous degree. 

Remark 1. Later on we will introduce additional variables A = (Ai, . . . , A„i) con- 
sidered as parameters, assign them appropriate weights and work in the extended 
ring C[a;, A] = C[xi, . . . , Xn, Ai, . . . , A^]. Even in the symmetric case the weights of 
the parameters will in general be different from 1. 

The Euler field associated with the weights Wi, . . . ,Wn is the derivation X = 
J2'^i Xid/dxi of C[x]. By construction, Xf = rf, r = deg/ G Q, for any quasiho- 
mogeneous polynomial / (the Euler identity). 

We put deg dxi = deg Xi = Wi . This extends the quasihomogeneous grading on 
all fc-forms: in the symmetric case, the degree of a polynomial /c-form will be k 
plus the maximal degree of its coefficients. Obviously, deg w = deg djuj for any form, 
provided that duj ^ Q. The Lie derivative Xlj of a quasihomogeneous form uj of 
degree r by the Euler identity is ruj. Note that degw > for all A:-forms with fc ^ 1. 

The norm of a polynomial in one or several variables is defined as the sum of 
absolute values of its (real or complex) coefficients. This norm is multiplicative. 
The norm of a A;-form by definition is the sum of the norms of its polynomial 
coefficients; it satisfies the inequality ||w A r]\\ ^ ||a;|| • \\r]\\ for any two forms w, r}. 

The exterior derivative operator is bounded in the sense of this norm if the degree 
is restricted: \\dbj\\ ^ (maxiWj)degw • In particular, in the symmetric case 

||cL;|| < ''ll'^ll, r = degw. Conversely, a primitive of an n-form /u can be always 
chosen bounded by the same norm 

Unless explicitly stated differently, a monomial (monomial form, etc) has always 
the unit coefficient. 
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3.2. Parameters. We will systematically treat the case when all objects (forms, 

functions etc.) depend polynomially on finitely many additional parameters A = 
(Ai, . . . , Am). Wo will denote by yl'^[A], fc = 0, . . . , n, the collection of A:-forms whose 
coefficients polynomially depend on A. For instance, the notation r] e A"~^[A] 
means that ry = ai{x, X) dxi A - ■ ■ A dxi A • • • A dXn with polynomial coefficients 

tti e C[a;, A]. 

In such case the norm of forms, functions etc. will be always considered relative 
to the ring C[x, A], that is, as the sum J2i absolute values of coefficients 

of the complete expansion in x, A. If the parameters A^ are assigned weights, we 
take them into account when defining the degree of the form. To stress the fact 
that the norm is computed relative to the ring C[.t, A] and not to C[a:] (i.e., that 
the situation is parametric), we will sometimes denote the norm by || • ||a. For an 
instance, ||2Aia;i|| = 2|Ai| ^ 2 = ||2Aia;i||A. 

3.3. Division by a quasihomogeneous difFerential df. The division modu- 
lus. If / e C[a;i,... quasihomogeneous polynomial having an isolated sin- 
gularity at the origin, then the multiplicity / of this singularity can be easily found 
by Bezout theorem, since no roots of the system of algebraic equations df /dxi — 0, 
i = l,...,n, can escape to infinity. In the symmetric case / = (deg/ — 1)". 
Choose any monomial basis (fi,. . . of the local algebra C[[xi, . . . , Xn]]/ (df), 
(df) = {-§^, ■ ■ ■ , ■§^)- Then the monomial n-forms /Xj = dxi A - ■ ■ A dXn form a 
basis of = A"/df A A"~^ over C: any n-form fj, can be divided out as 

I 

IJ. = ^Cil^i + df Ar], CiSC, T?eyl""\ (3.1) 

i=l 

with appropriate constants Ci, . . . , c/ S C (coefficients of the "remainder" ^ Cj/Uj) 
and a polynomial form rj S A"~^ (the "incomplete ratio"). Moreover, if /i is 
quasihomogeneous, then the decomposition (3.1) contains only terms with deg/itj = 
deg fj, and degrj = deg/x — deg/. This immediately follows from quasihomogeneity 
and the uniqueness of the coefficients q. Form this observation we also conclude 
that all monomial forms of degree < deg / must be among /Xj, and, moreover, any 
monomial form of degree greater than maxjdeg/Xj, is divisible without remainder 
hjdf. 

The choice of the monomial forms /Xi spanning the quotient, is not unique, though 
the distribution of their degrees is. Denote by p = p(/) the maximal difference 

p{f) = max deg /x, — min deg /Xj = max deg (pi — min deg ipi . (3.2) 

i i i i 

The following results are well-known. 

Proposition 1. 1. In the symmetric case p{f) < I — {r — 1)" [AGV85, §5.5]. 

2. In the bivariate case n = 2 the inequality p{f ) < r = deg / holds if and only 
if f is a "simple singularity" of one of the following types, 
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see e.g., [AGV85, §13, Theorem 2]. □ 

Prom these observations it can be immediately seen that the division with re- 
mainder (3.1) is a bounded linear operation in the space of all n-forms of restricted 

degrees. 

Lemma 1. Assume that f & is a quasihomogeneous polynomial having an 
isolated critical point of multiplicity I at the origin, and the monomial n-forms 
/Xi, . . . , /ij G A", form the basis of yl". 

Then there exists a finite constant M < +oo depending only on f and the choice 
of the basis {^i}, such that any n-form /j. G A" can be divided with remainder by 
df as in (3.1) subject to the following constraints, 

degr?<degM-deg/, ||7?|| + ^ |ciK M||m||. (3.3) 

// the form, fi is quasihomogeneous, then deg r] = deg /U — deg / and Ci can be nonzero 

only i/dcg/A; = dcg/i. 

The constant M depends on the choice of the monomial basis {i^i}. The optimal 
choice of such basis (out of finitely many possibilities) results in the smallest value 
M = M{f) that depends only on /. We will always assume that the basis {/ij} is 

chosen optimal in this sense. 

Definition 1. The minimal constant M{f) corresponding to an optimal choice 
of the monomial basis of the quotient AJ, is called the division modulus of the 
quasihomogeneous polynomial / e A*^. 

Corollary 1. Assume that /i G A" [A] depends polynomially on additional parame- 
ters A. Then fi can be divided with remainder by df so that the remainder and the 
incomplete ratio depend polynomially on A with the same division modulus, 

Ci = Ci{X) e C[\], i = l,...,n, r?GA"-^[A], 

IH| + ^lhHM(/)|H|, ||-|| = ||-IU. 

Proof of the corollary. Every monomial from the expansion of fj. in x,X can be 
divided out separately by df which is independent of A. □ 

Proof of the Lemma. Let M be the best constant such that (3.3) holds for all mono- 
mial n-forms with deg/x ^ L It is finite since there are only finitely many such 
forms. In particular, since any form of degree I is divisible by df by Proposition 1, 
the respective fraction 77 will be of the norm at most Af||/i||. 

Writing an arbitrary monomial n-form of degree > I as a product of a monomial 
form of degree I times a monic monomial function x" G C[a;], a G Z", we con- 
struct the explicit division formulas (without remainders) for all monomial forms 
of higher degrees. The division constant will be given by the same number M, since 
multiplication by a monic monomial preserves the norms of both \\fi\\ and ||?7||. 

All the other assertions of the Lemma are well-known [AGV85]. □ 

3.4. Computability of the division modulus. Despite its general nature, the 
above proof is constructive, at least in the low dimensional cases n = 1, 2, allowing 
for an explicit computation of the division modulus in these cases. 

The one-dimensional case is trivial: for the monomial f{x) = x^ the division 
modulus M(/) is equal to r and it can be obviously recalculated for any other 
principal homogeneous part. The "special case" of a multivariate polynomial f{x) = 
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Xi-\ see [GIOl], is reducible to the one-dimensional situation. In this case 

I = (r — 1)" monomial forms x" dxi A • • • A dxn with ^ ^ r — 1 form the 
basis, and the corresponding division modulus is again equal to r. This example 
admits an obvious generalization for quasihomogeneous "special polynomials" with 
different weights. 

For a bivariate truly homogeneous polynomial / (i.e., in the symmetric case, 
the most important for applications), the division modulus M for all higher de- 
gree forms (dcg /i ^ 2 dcg /) can be explicitly computed as the norm of the inverse 
Sylvester matrix for the partial derivatives and [NYOl]. The " quasimonic" 
polynomials, introduced in that paper, are defined by the condition M(/) = 1, 
which in many respects is a natural normalizing condition for multivariate polyno- 
mials. 

The choice of the basic forms even in the symmetric bivariate case depends 
on /: while it is gcncrically possible to choose tlicni as x"^x'^'^ dxi A dx2 with 
^ Q;i,2 ^ J' — 1, for a badly chosen / some of these forms of degree greater than 
r = degf can become linear dependent modulo df, requiring a different choice. In 
order to avoid making this choice, one may allow a redundant (i.e., linear dependent) 
collection of generating forms /Xj. Choosing all monomial forms of degree < 2r 
makes the corresponding division for low degree forms trivial, so that the division 
modulus M{f) is determined only by division of forms of higher degree. Details 
and accurate estimates in the bivariate symmetric case can be found in [NYOl]. 

To describe the division modulus M(/) in the case of n ^ 3 variables is a 
considerably more difficult problem, though it still can be reduced to analysis of 
finitely many monomial divisions. One can (at least, theoretically) express M{f) 
via lower bounds for minors of certain explicitly formed matrices. 

Remark 2. It is worth mentioning that the division modulus M(/) is not directly 
related to the norm ||/||, even in the symmetric bivariate case. If deg /i ^ I and /i = 
df A r?, then ^ \\df\\ \\ri\\. On the other hand, ^ M-i||r?|| by the definition 
of M(/). Therefore M(/) ^ \\df\\-^ > ||/||-S that is, the division modulus for 
a polynomial / with the small norm must be large. The inverse is not true: a 
polynomial with a small division modulus can have a very large norm. Simple 
examples can be constructed in the form f{x) = c'[l^{xi — XiX2) with sufficiently 
close values of the parameters Aj e [0, 1] and a suitably chosen normalizing constant 



3.5. Division by /. We begin by establishing an analog of the Euler identity 
in the Brieskorn module. It plays the central role for explicitly constructing the 
decomposition (2.6). 

Lemma 2. Assume that f £ A" is a quasihomogeneous polynomial of degree r. 
Then any polynomial n-form divisible by df in , can itself be divided by f in the 
Brieskorn module Bf. It also admits a polynomial primitive divisible by f. 

In other words, for any form r/ € there exist four forms fj, G A", u e 

and ^, ^' e such that 



The degrees of all forms are all equal to degrj in case the latter is 

quasihomogeneous. 



ceC. 



df Ar, = ffi + df AdC 

= d{fuj) + df Ad^'. 



(3.4) 
(3.5) 
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The division operation is always well-posed in the sense that the decomposition 
(3.5) can be always chosen to meet the inequality 

M + U'W ^ (n + 3)degr7.||r7|| (3.6) 

(a similar inequality can be proved also for the first decomposition (3.4)). 

Proof. Note that for any n-form fi G A" and any vector field X on C", 

{X.f)l^ = (ix df) M = d/ A ^, 

where is the inner antiderivative, since d/ A /x = 0. We wiU need this formula for 
the case when X is the Euler vector field. 

To prove the first divisibility assertion (3.4), we have to show that the identity 

df Arj = fn + df Ad^ (3.7) 

can be always resolved as a linear equation with respect to fx and ^ for any choice 
of r]. Using the Euler identity for functions and the above remark, we represent ffi 
as a form divisible by df, 

/m = r-'{Xf) n = r-\i, df)fi = dfA r'^ i, (3-8) 

The equation (3.7) will obviously be satisfied if 

that is, when ry is cohomologous to i^ ji. This last condition is equivalent to the 
equality between the exterior derivatives 

dr] = r~^ di^ 11 = r~^Xn, 

since by the homotopy formula, di^ ^ = X^i — i^d^x = Xfi. Thus resolving the 
equation (3.7) is reduced to inverting the Lie derivative X on the linear space of 
n-forms. 

We claim that the linear map /i Xfi of A" to itself, is surjective (and obviously 
degree-preserving) , guaranteeing thus solvability of the last equation for any choice 
of 77. Indeed, any monomial n-form /i^ = x" dxi A • • • A dxn is an eigenvector of 
X with the strictly positive eigenvalue deg Hc, n (recall that the weights Wi are 
normalized so that the volume form dxi A ■ ■ ■ A dxn is of degree n). Thus X is 
surjective on yl" (actually, bijective) and one can choose /x = rX^^{dr]). The norm 
of the inverse operator X~^ does not exceed (r/n) degt] in the symmetric case. The 
proof of (3.4) is complete. 

To prove the second assertion (3.5), we transform it using (3.8) as follows, 

df Ar] = fduj + df A{uj + d^) = r-^df Ai^ dco + df A (lu + d^), 

which will be obviously satisfied if 

r) = r~'^ d(j + (J + d^ . (3.9) 

Taking the exterior derivative as before, we reduce this equation to the form 

dr] = r~^ di-ji du + du = r~^ Xjj, + 11, fi = du. 

Solvability of this equation with respect to /i (and hence to uj) for any left hand 
side dr] follows from invertibility of the differential operator r~^X + 1 on the linear 
space of polynomial n-forms (1 stands for the identity operator). Exactly as in 
the previous situation, all monomial n-forms are eigenvectors for {r~^X + l)|yin 
with the positive eigenvalues, all greater or equal to r~^n -|- 1, hence r~^X -|- 1 is 
invertible on and oj can be chosen as a primitive of {r~^X + l)~^dr]. 
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To prove the inequality between the norms, notice that fj, = duj satisfies the 
inequality ^ \\dri\\ ^ dcg?] ||77||. A primitive u! can be always take of the norm 
llwll ^ \\duj\\. Together this yields < degry ||r;||. 

The norm ||^'|| can be found from (3.9). Clearly, || i^ fj,\\ < n||/u|| because of the 
choice of the weights dcg Xi which satisfy the condition = n. Substituting 

this inequality into (3.9), wc obtain 

lle'll < IM^ II ^ ll'7ll + n\\M\ + M < (n + 2) degT? ||7?||, 
since degw = degry > 1. □ 

3.6. Generating Petrov and Brieskorn modules: the algorithm. Division 
by the gradient ideal together with the Euler identity as formulated in Lemma 2, 
allows for a constructive proof of the representation (2.6) for an arbitrary semi- 
quasihomogcncous polynomial F. 

Let F = f + h £ C[xi, . . . , Xn] be a semiquasihomogeneous polynomial with the 
principal quasihomogeneous part / and the lower-degree part h. Denote as before 
by /ii, . . . , ^/ e A" the forms spanning — A"/df A A^~^ (note that the quotient 
is computed using only the principal part /). We claim that: 

1. any n-form /x G A" can be represented as 

; 

M = ^aiMi + «^^AdC, G C[F], C e (3.10) 

i=l 

2. any (n — l)-form u) G A^~^ can be represented as 

I 

w = Y,Pi^i+dF A^ + d^, ft gC[F], e,C'e^""^- (3.11) 

The construction of the decomposition (3.10) begins by division of /x by df as 
explained in Lemma 1: 

H = ^Ciiii + df AT], CjGC, r/GA"-^ 

If deg/i < r = degf = degF, then the incomplete ratio is in fact absent, r] = 0, 

and wc arrive to a particular case of (3.10) with qi = Ci of degree (constants). 

If deg/x is higher than r, we transform the term df Arj using Lemma 2 and then 
substitute f = F — h: 

M - Czl^i = /m' + df AdC 

= Fn' + dF AdC- fi", fi" = hfi' + dhA d^. 

Obviously, both fi' and jj," are of degree strictly inferior to deg/x, which allows 
to continue the process inductively. Assuming that the representations (3.10) are 
known for both /i' and /i", wc substitute them into the last identity and after 
collecting terms arrive to a representation for /x. In the symmetric case the inductive 
process cannot take more than deg /x — r steps. It is a direct analog of the process 
of division of univariate polynomials, see also [NYOl]. 

To construct (3.11), we divide dtu by df. If degw < r, then the incomplete ratio 
is absent and we obtain a special kind of (3.11) exactly as before. 
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Otherwise in the division with remainder 

du = ^Ciduii + df A T], Ci e C[A], rj e A"~^[A], 

substitute df Ar] = d{fuj') + df Ad^ and pass to the primitives. We obtain 

a; - ^ CiUJi = fJ + df A^ + d^' 

= Flo' + dF A^ + d^' - J', J' = hu' + dhAC 

For the same reasons as before, the degrees of w', oj" are strictly smaller than degw, 
hence the process can be continued inductively. 

Remark 3. In a somewhat surprising way, it turned out impossible to transform 
directly the decomposition (3.10) for the form du S A" into (3.11) for w. 

3.7. Effective decomposition in the Petrov module. The construction above 
is so transparent that any qualitative as well as quantitative assertion concerning 
these expansions, can be immediately verified. 
Wc will show that 

1. all terms of the decomposition (3.11) depend polynomially on the lower order 
terms of F, assuming that the principal part if fixed, and 

2. the wcU-poscdncss of the construction is determined solely by the division 
modulus M{f) of the principal homogeneous part. 

In order to formulate the result, consider a general semiquasihomogeneous poly- 
nomial with the prescribed principal quasihomogeneous part, 

F{x,X)=f{x) + h{x,X), h{x,X)= ^Js{x), (3.13) 

deg/<,<deg/ 

where fi, . ■ ■ ,fm € C[a;i, . . . , Xn] are all (monic) monomials of degree strictly infe- 
rior to r = deg/, arbitrarily ordered. We treat the coefficients Ai, . . . , as the 
parameters of the problem, assigning to them the weights so that 

deg As + deg fg = deg / = r for all s. 

This choice makes the entire polynomial F quasihomogeneous of the same degree 
r in the ring C[a;, A] = C[a;i, . . . , a;„, Ai, . . . , A^]. Instead of the ring C[F], the 
coefficients pi of the decomposition (3.11) will belong to the ring C[F, X] and their 
quasihomogeneity will be understood in the sense that the formal variable F is 
assigned the weight deg F = r. 

Theorem 1. If the quasihomogeneous polynomial f € C[x] has an isolated critical 
point at the origin and F e C[x, X] is a general semiquasihomogeneous polynomial 
(3.13), then any polynomial quasihomogeneous {n — l)-form to e A"~^[A] of degree 
k can be represented as 

I 

uj = YpiUJi + dF A$,+d£,'. (3.14) 

The coefficients Pi £ C[-F, A] and the (n — 2)-form,s G yl"~^[A] are all polyno- 
mial and quasihomogeneous jointly in F, X {resp., in x, A) of the degrees k — degWj, 
k — r and k respectively. 
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The norm of the coefficients relative to the ring C[F, Ai, . . . , A^] is explicitly 
bounded in terms of n,r,k and the division modulus M{f). In particular, for the 
symmetric case when degxi = • • • = dega;„ = 1, 

I 

WViW ^ k\ r'=("+3)M'= \\uj\\, k = degu;, M = M{f), || • || = || • lU- (3.15) 

i=l 

Remark 4. The fact that the form to is quasihomogeneous, is not important: any 
polynomial form is the sum of quasihomogeneous parts, each of them being divisible 
separately. 

Remark 5. Even in the symmetric case, the degrees of the parameters are different 
from 1: deg Xg = r — deg fg will take all natural values from 1 to r. 

Proof of the Theorem. The first assertion of the Theorem (on polynomiality and 
quasihomogeneity) follows from direct inspection of the algorithm described above, 
since all transformations on each inductive step (exterior differentiation, division 
by df which is independent of A, and the Euler identity in P/) respect the quasi- 
homogeneous grading. 

The only assertion that has to be proved is that on the norms. In order for a 
sequence of increasing with k real constants Cfe > to be upper bounds for the 
decomposition (3.14), 

; 

\\pi\\ < Ck ||w||, for all oj with deg a; < k, 

i=l 

they should satisfy a certain recurrent inequality which we will instantly derive 
from the suggested algorithm. 

Denote by pi E C[F, A] (resp., by p'^ and p") the polynomial coefficients of the 
decomposition of the forms ui (resp., ui' and iv") from the identity (3.12): since the 
degrees of both u}',ui" are less than k and the sequence Ck is increasing, we have 

J2\\p'i\\^Ck-i\\io'\\, 53iiKIKCfe_i|K||. 

i i 

Multiplication by F corresponds to a shift of coefficients in the decomposition of 
w'. Thus from (3.12) follows the inequality 

E \M ^ E ii^^ii + E \\Pi\\ + E ii^'"!! < E ii^^ii + ^'^-idi'^'ii + 

i i i i i 

By Lemma 2, ||a;'|| < (n + 3)A; ||?7||. The norm of the inferior part h is by definition 

equal to the number of terms, that is, the number of monomials in n variables of 
degree < r — 1. Therefore \\h\\ ^ r" and \\dh\\ ^ r"+^. This implies an upper bound 
for ||u;"||: 

< ll^ll + \\dh\\ m\ ^ m + WdhmWW + \m < 2r"+^(n + 3)^ M 

by Lemma 2. Finally, \\ri\\ +Y1 Ik-ili ^ ll'^ll by definition of the division modulus 
M = M{f). Assembling all these bounds together, we conclude that 

E IIPilK Arilwll + Cfc_i • 3r"+i(n + 3)fc ||a;||. 
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Thus the increasing sequence Cfe ^ 1 will form upper bounds for the norms of the 
coefficients of decomposition for polynomial forms of degree < k, provided that 

Cfc > AkCk-i, A > 4r"+^(n + 3)M > r"+^M 

(notice that r ^ 2), which can be immediately satisfied if we put 

This proves the inequality for the norms. □ 

Note that the bound established in this Theorem, is polynomial in M = M(/) 
and (for a fixed r) factorial in A; = degu), that is, only slightly overtaking the 
exponential growth. 

3.8. Nonhomogeneous division. By a completely similar procedure one can de- 
scribe the result of division by a nonhomogeneous differential dF as a sequence of 
divisions by the principal homogeneous part df. 

More precisely, if e ^"[A] is a polynomial n-form polynomially depending on 
the parameters Ai, . . . , Am and F = f + '^sfs is as in (3.13), then there exists a 
representation 

IJ, = Y,Ci{X)lJ^i + dF Ar], ci,...,c„ e C[A], r? e ^"-^[A], (3.16) 

polynomially depending on parameters. If /i is quasihomogcncous, then so arc Ci 
and r], with deg Cj = deg /x — deg jJi and deg rj = deg /x — deg F. Moreover, the ratio 
(IIciIIa + II'7||a)/||a*||a is bounded in terms of deg/i and the division modulus M(/) 
of / only. 

Indeed, dividing jj by df yields 

A« = ^ CiHi + df At] = CiHi + dF Arj- iJ,', fi' = dh A r], 

where h = F — f, hence deg/i < deg/ = degF and therefore degfi' < dcgii. 
This means that the process of division can be continued inductively. Since ||/i'|| ^ 
\\h\\ \\r]\\ ^ constr,n M{f) the norms of the remainder and the incomplete ratio 
are bounded in terms of M{f) and the degrees. In the symmetric case the bound 
looks especially simple. 

Proposition 2. In the symmetric case of all weights equal to 1, the division of a 
form of degree k = deg is bounded as follows, 

I 

Mx + E ll'^^IU ^ Mk{F) ■ M, M,{F) = fcr"('=-'') {M{f))K 

Proof In this case \\h\\ ^ r", so that ||^'|| ^ Mr"||^||, and finally ||77|| + J2 < 
M||/il|(l + K + ■■■ + Kdeg/x-r-)^ ^Yieie K = Mr"". Thus the norm of the non- 
homogeneous division operator obviously does not exceed M*'(A;r"^'°~''^). This ex- 
pression is exponential in fc = deg/i and polynomial in M = M{f). □ 
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4. PiCARD-FUCHS SYSTEM FOR AbELIAN INTEGRALS 

Consider a quasihomogeneous polynomial / e A° of degree r = deg/ having 

an isolated singularity of multiplicity I at the origin. As before, let fii,...,iii 
be generators of over C and u)i,...,LUi their monomial primitives. Consider 
the general semiquasihomogeneous polynomial F = f + J2T "^^fs S C[a;, A] as 
in (3.13) with the fixed principal part /, whose coefficients Ai,...,Am arc the 
natural parameters. Consider in the parameter space the locus S such that 
for A e C™ \ S the level set {x € C" : F{x, A) = 0} is a nonsingular algebraic 
hypersurface. Denote by -T = r{X), X <^ E, any continuous family of (n — l)-cycles 
on the zero level. The Abelian integrals 

= f CO,, i = l...,l (4.1) 

Jr{\) 

are well defined multivalued analytic functions on C" \ S. In this section we will 
derive a Pfaffian system of linear equations satisfied by these integrals. 

We will always assume that the weights of the parameters Xg are chosen so that 
F becomes a quasihomogeneous polynomial in x, X of degree r: deg Xg — r — deg fs- 
The enumeration of the monomials fs begins with the free term /i = 1 of degree 
so that the respective coefficient Ai is necessarily of degree r. Recall that p{f) is 
the maximal difference (3.2) between the degrees of the forms /i^. 

Theorem 2. There exist {I x l)-matrix polynomials Co(A),Ci(A), 

Co(A) = Ai •1 + C'(A2,...,A„), 

deg Co ^ r + p(/) , deg ^ deg /, + p(/) , s = 1 , 

{the degrees are quasihomogeneous) , such that on C™ \ S 

^ (Co(A)J) = Cs{X)I, s = 1, . . . , m. (4.3) 

The norms ||Cs||a d^e bounded by a power of the division modulus M{f). 

In other words, the column vector function 7(A) on the complement to S satisfies 
the matrix Pfaflaan equation 

dl = 07, O = Cq"^ ■ (-dCo + ^Cs dX^ , (4.4) 

^ s=l ^ 

with a rational matrix-valued 1-form O having the poles only on the locus Z" = 
{detCo = 0} C C™. Here d is the exterior derivation with respect to the variables 
A. only: for c(A) £ C[A], dc = ^ 

The proof is constructive. The description of the matrix polynomials Cs{X) is 
given below. 

4.1. Gelfand— Leray derivative with respect to parameters. 

Lemma 3. If to G A"~^ is a polynomial form with constant {independent of A) 
coefficients, and rjs G yl"~^[A] any form satisfying the identity 

fsdoj = -dF Arjs, (4.5) 



, ^^(A), 



(4.2) 
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{recall that fs = ■§&), then 



O'^s Jr{\) Jr{\) 

Proof. To derive this formal identity, we express = H{x) from the equation 
F{x,\s) = 0, assuming all other parameters fixed, and apply the Gelfand-Leray 
formula to H: for (4.5) to hold, it would be sufficient iirj = rjs satisfies 

duj = dH A Tj. 

It remains to observe that by the implicit function theorem and the definition of 
the parameters, 

dF dF 
dF+^dH = 0, ^=fs. 

OAs OAs 

Here and above d stands for the exterior derivative with respect to the "spatial" 
variables xi,. . . ,Xn- □ 

The standard Gclfand Lcray derivative appears for the parameter occurring be- 
fore the constant term fi = 1 (modulo the sign). 

4.2. Derivation of the system: beginning of the proof of Theorem 2. Di- 
vide each of the forms FfXi G -^"[A], /Xj = dcot, by dF with with the remainder 
coefficients and the incomplete ratios polynomially depending on A as in Proposi- 
tion 2: 

I 

Fjjti = dF Arii + 'Y^ Cij fij , Cij = Cij (A) . (4.6) 

i=i 

Clearly, the quasihomogeneous degree deg Cij in C[A] is equal to r+deg /Xi — deg rrij ^ 
p{f) + ^ (cij = if the difference is negative). 

Let Co = Co(A) be the (Z x Z)-matrix polynomial with the entries Cij{X). Since 
dF does not depend on Ai (the free term of F), while the only term depending on 
Ai in F^i is XiHi, the dependence of Co on Ai can be immediately described: the 
corresponding remainder coefficients (Ai) for the division of Ai/Xj by dF form the 
scalar matrix Ai • 1 (the incomplete ratio is absent). 

Since Cij do not depend on x (being "constants depending on the parameters"), 
the identity (4.6) implies that 

d{Fuii - CijOJj) = -dF A (-Wj - r]i), i = l,...,l. 

Let 

i^is = -fsii^i+Vi), i = l,...,l, s = l,...,m. 

All these forms are polynomial and polynomially depending on parameters. Their 
degrees can be easily computed: deg rji = deg m = deg Wi, deg w- = deg /s+deg /i^. 

By the parametric Gelfand-Leray formula (Lemma 3), the partial derivatives of 
integrals of the forms Fwi — Cijivj over the cycle r{X) C {F = 0} C C" are 
equal to the integrals of the forms Since the terms FuJi vanish on r{X) for all 
values of A, we have 



3 
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The forms were chosen to generate the Petrov module Pp over C[F, A], so each 
of the AbeUan integrals § i^l g ^e expressed as a polynomial combination, 

I 

for all i, s. Denote by Cg = Cs{X) the polynomial {I x Z)-matrix function formed by 
the free terms of the polynomials Pij,s{'j X)' 

C«(A)= [PijAF,>^)\F=o\[j=i' s = l,...,m. 

All other terms, being divisible by F, disappear after integration over the cycle 
on the level surface {F = 0}. Collecting the terms, we conclude that the partial 
derivatives of the column vector function J(A) = (/i(A), ...,/;( A)), = ^ uii, we 
have 

^ = CJ, s=l,...,m. 

4.3. Bounds for the norms: end of the proof of Theorem 2. The con- 
struction described above, does not yet imply the assertion on the norms of the 
matrix polynomials Co, ... , Cm for only one reason: multiplication hy F = f + h, 
h = X^As/s, is not a bounded operator. While multiplication by h increases the 

norm at most by \\h\\\ = constn,,. (not exceeding (r — 1)" in the symmetric case), 
the norm ||/|| cannot be bounded in terms of M{f), as required in the Theorem 
(see Remark 2). 

To correct this drawback, exactly as in [NYOl], the division line (4.6) should be 
first prepared using (3.8) as follows, 

Fidi = {f + h)fi^=dfAr]- + hfi, = dFAf]l+ii'i, 

(4.7) 

Vt^r~'^i^^i, H', = hfit - dh At]'^, 

where (we again make all estimates for the symmetric case only), 

IKIK Wr)||Mi||, \\^^[\\^\\h\\{l + r){n/r)\\^,,\\. 

Then forms /i'^ should be divided by dF with remainder: since their norms are 
bounded by a constant depending only on n, r (the norms of the monomial forms 

are equal to 1), the results of such division will be bounded by suitable powers 
of M{f) by virtue of Proposition 2. 

Collecting the terms, we conclude that the coefficients Cij S C[A] of the cor- 
responding remainders in (4.6) and the incomplete ratios rji € A"~^[A] will be 
bounded by expressions polynomial in M{f). 

The rest of the derivation remains unchanged and the estimates completely 
straightforward: the polynomial bounds for rji imply those of the polynomial coeffi- 
cients pij^s G A] by Theorem 1. This proves the last assertion of Theorem 2. □ 

5. Observations. Discussion 

The algorithm of derivation of the Picard-Fuchs system in the PfafHan form is 
so transparent that many things become obvious. 
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5.1. Bounds. Though the matrix polynomials Cs{\) are not quasihomogeneous 

(their entries have different degrees), the determinant detCo(A) is a quasihomoge- 
neous polynomial from C[A]. Its degree can be immediately computed as Ir from 
the explicit representation (4.2). This same representation proves that this deter- 
minant, equal to A" + polynomial in(A2, . . . , A^), does not vanish identically, so 
that the system (4.4) is indeed meromorphic. 

Moreover, the norm of the inverse matrix Cq^ can be explicitly majorized in 
terms of the distance to the critical locus. One possibility to do this is to con- 
sider the sections Ai = 1 and apply the Cartan inequality as in [NYOl], using the 
quasihomogeneity. 

5.2. Spectrum. The spectrum of Co(A) can be also easily computed: it consists 
of all I critical values of the polynomial F{x,X), at least when F{-,X) is a Morse 
polynomial. To see this, it is sufficient to evaluate both parts of (4.6) at any of 
I critical points ai,...,ai G C". The column vectors Vi = {(fi{ai),...,(fi{ai))'^, 
i = 1, . . . , Z, are the corresponding eigenvectors (recall that /x, = ipi dx\ A • • • A dx„). 

5.3. Hypergeometric form. Restricting the PfafSan system (4.4) on the one- 
dimensional complex lines Ag = const, s = 2, . . . , m, parameterized by the value of 

t = Ai, one obtains a parameterized family of Picard-Fuchs systems of ordinary 
differential equations. In this case only the matrix C\ is relevant. 

By Theorem 2, it is quasihomogeneous of degree < p{f) jointly in the variables 
Ai,...,A„i. If p{f) < r = degAi, then Ci cannot depend on Ai and hence the 
Picard-Fuchs system in this case will have the hypergeometric form (1.1). By 
Proposition 1, this happens only when / is a simple quasihomogeneous polynomial 
of one of the types listed there. For hyperelliptic polynomials (the singularity of 
the type Ak) this was well-known, see [NYOl]. In turn, the hypergeometric form 
implies that all singular points of the Picard-Fuchs system are Fuchsian. 

5.4. Logarithmic poles. For the full PfafBan system (4.4) the polar locus, occur- 
ring where detCo(A) vanishes, is of multiplicity 1 (it is sufficient to produce just 
one value of the parameters A such that i^(-,A) has simple critical points). Yet it 
is not the characteristic property. 

A rational 1-form ui analytic outside a hypersurface S' = {g = 0} c C", g being 
a polynomial without multiple factors, is said to have a logarithmic singularity on 
this hypersurface, if both gu} and dg Alu extend as polynomial forms across S' on 
C™. 

This is only one of several close but non-equivalent definitions, probably the 
strongest possible. It ensures that the restriction of u on any holomorphic curve 
7 cutting S' at a point a, has a Fuchsian singularity with the residue independent 
on the choice of 7, depending only on the point a. 

The basic question concerning the system (4.4) is whether this system itself or a 
suitable gauge transformation of this system with a rational matrix gauge function, 
are Fuchsian with bounded residues. If the answer is positive, this would mean a 
positive solution of the infinitesimal Hilbert problem. 

Using symbolic computation for implementing the algorithm, we discovered that 
in the hyperelliptic case (singularity of the type Ak) the Picard-Fuchs system (4.4) 
indeed has only logarithmic poles until the degree fc = 6 of the polynomial / = 
Xi + This naturally suggests the following conjecture. 
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Conjecture. All singularities of the Picard-Fuchs system (4.4) are only logarith- 
mic poles on E' = {det Co = 0}. 

It would be interesting to verify this conjecture for other simple singularities 
listed in Proposition 1, perhaps first by symbolic computation. 

The next step could be to study the behavior of residue of (4.4), the matrix 
function defined on the regular part of E' , checking whether it is bounded near 
singular points of the discriminant. 

5.5. Singular perturbations. The polynomial dependence of the matrices Cg on 
the lower degree coefficients of the polynomial F = f + ■ ■ ■ fails for the coefficients 
of the principal part. Though apparently rational, this dependence certainly must 
exhibit singularities when / degenerates into a quasihomogeneous form with non- 
isolated singularities. The Picard-Fuchs system in such cases may have singular 
points corresponding to atypical values of F. Their appearance must somehow be 
related to the fact that the division modulus explodes when such degeneracy occurs, 
thus creating a singularly perturbed system of linear differential equations. These 
phenomena seem to be worth of detailed study. 
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